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-  >  To  provide  a  simpler  procedure  for  performing  the  required  integration 
and  thus  make  the  present  method  more  rapid  and  convenient,  the  virtual 
height  data  have  been  modeled  in  terms  of  a  polynomial  governed  by  a  single 
parameter.  The  parameter  was  determined  by  comparison  of  the  model 
with  the  data.  The  results  using  the  model  were  in  good  agreement  with  those 
obtained  by  direct  use  of  the  data. 

v  The  present  method  affords  a  simple  and  convenient  means  of  obtaining 
real  height  from  virtual  height  data,  within  the  accuracy  of  the  experimental 
data. 
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Application  of  a  Direct  Method  for 
Real  Height  Determination  to  Two  lonograms 


I  I M  Konl  <  TION 

in  a  previous  report,  1  a  direct  procedure  was  developed  for  the  determination 
of  real  neight  from  virtual  height  data  obtained  in  the  presence  of  a  magnetic  field, 
i'hc  method  was  based  upon  the  representation  of  the  dispersion  relation  between 
index  of  refraction  and  electron  density  py  a  simple  power  law,  resulting  in  an 
Abel  integral  equation  for  the  relation  between  real  height  and  virtual  height,  which 
could  then  be  solved  explicitly  for  the  real  height  as  a  function  of  frequency.  To 
study  the  accuracy  of  the  method,  two  parabolic  models  of  electron  density  were 
convened  into  virtual  height  data,  lrom  which  the  real  height  was  determined  by 
the  present  method.  The  results  obtained  were  in  good  agreement  with  the  original 
models . 

In  the  present  report,  the  method  has  been  extended  to  the  analysis  of  field 
data.  To  assess  the  accuracy  and  reliability  of  the  method,  ionograms  have  been 
chosen  for  which  accurate  numerical  calculations  of  real  height  were  available  for 
comparison. 
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Boi- *  dec ,  dorado,  '  and  Ga;  ir  ,  ?i  anct  ‘‘  aia  -iiow.i  ir  Figures  1  and  2. 

the  ;0'.v<-Tjt  i.-.-quenoy  1  .  tot  *  >»:  'h  iwsond*  u.. ;«•  is  sva .  -able  is  usually  limited 

to  a  ranee  c:  1  to  ?  MHz,  it  is  no  cossa:—  to  clete  -Mime  the  initial  height  h ,,  below 

which  the  etfect  of  ionization  mav  be  ccnsirterc  .  negligible.  To  determine  it 

is  necessary  to  assume  a  variation  in  virtual  hem-u  h‘  iron  f  down  to  the  fre- 
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quenc'.  f  corresponding  to  the  initial  height.  Ovtr  c  the  Gar env  ionogram  is  llatter 
in  the  lower  frequency  region  than  the  Boulder  curve,  we  may  anticipate  that  the 
initial  height  calculation  will  be  more  accurate  for  Garchy  than  for  Boulder. 

As  shown  in  Reference  1,  the  real  height  is  retained  from  the  integral 
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where  a  is  a  frequency  dependent  parameter, 
sin  (ir^-  ) 
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and  fv  is  the  maximum  probing  frequency.  Since  the  lowest  frequency  for  which 
data  is  available  is  f  ,  we  write  Eq.  (l)  in  the  form 

f 

r  1  h’Cf>  -  h0 

h  ‘  h0  =  J  pfo)  fl-o  (f2  _f2f/ 2  df  • 2  (3) 

f0  V 

so  that,  the  integral  gives  the  increment  in  real  height  measured  from  h0-  A  curve 
of  the  parameter  a  as  a  function  of  frequency,  taken  from  Reference  1,  is  given  in 

Figure  3. 
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Figure  l.  Boulder,  Colorado,  Ionogram, 
20  February  197  0,  1900:27  MST.  Declina¬ 
tion  I  =  67. 5° 
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Figure  2.  Garchy,  France,  Ionogram, 
10  March  1901,  02 00  Local  Time. 
Declination  1  =  03° 
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Figure  3.  Plot  of  Parameter  a  Against  Frequency;  6  is  Angle  Between 
Ray  and  Magnetic  Field 


3.  DETERMINATION  OF  INITIAL  HEIGHT 


We  have  found  the  one-start  procedure  of  Howe  and  McKinnis4  convenient  to 
use  and  easy  to  adapt  to  our  method.  We  write  Eq.  (3)  in  the  form 
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where  x  -  f/fy  is  a  normalized  frequency  variable.  We  now  assume  a  linear  varia' 
tion  of  h'  between  xQ  and  nearby  points  Xj  and  x9  at  which  h'  is  known  and  write 
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4.  Howe,  H.  H.  ,  and  Mclnnis,  D.  E.  (ipt>7)  Ionospheric  electron-density  profiles 
with  continuous  gradients  and  underlying  ionization  corrections.  II.  Formu 
lation  for  a  digital  computer,  Radio  Sci.  2 (New  S.-riesHNo.  10);113f>-ll58. 
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whore,  since  o  is  small,  we  have  taken  a  -  t)  and  p(o)  =  1.  Calculations  show  that 
chanties  in  tor  small  a  produce  only  slight  changes  in  the  integral.  To  deter¬ 
mine  h  ,  we  reed  an  additional  relation  among  the  real  heights  h^,  h^,  and  h^  in 
1-qs.  (  -> )  and  (d).  it  e  shall  make  the  reasonable  assumption  that  these  real  heights 
he  or.  a  straight  line  and  write 
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i-  re  111  these  equations  we  obtain 
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Trial  calculations  showed  the  value  of  was  slightly  sensitive  to  the  choice 
of  x|(;  a  variation  of  0.25  in  yielded  a  change  of  about  3  km  in  hg  for  Boulder 
and  less  than  one  for  Garchy.  Since  the  calculated  value  of  h(1  for  Boulder  was 
somewhat  low,  a  parabolic  scheme  was  tried  in  place  ot  the  preceding  linear 
mi  thod.  In  plan-  of  l.qs.  (  >>  and  (,;)  we  write 
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while  Eqs.  (8),  (9),  and  (10)  are  replaced  by 
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The  parabolic  fit  was  found  to  substantially  improve  the  initial  height  for  the 

Boulder  ioriograrn  while  leaving  its  value  for  the  Garchy  case  virtually  unchanged. 

We  have,  therefore,  utilized  the  parabolic  fit  in  the  determination  of  initial  height. 

In  the  absence  of  data,  the  determination  of  h^  by  a  method  such  as  the  present 

one,  which  relies  on  ordinary  ray  data  only,  would  ignore  any  sudden  changes  in 

virtual  height  below  the  minimum  frequency.  If,  however,  data  for  the  extraord- 
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inarv  ray  is  also  used,  sudden  changes  in  the  real  height  can  be  detected.  ’  For 
example,  tht  calculated  value  of  real  height  for  the  Boulder  case  in  Reference  2 
shows  a  sudden  rapid  drop  when  the  frequency  decreases  to  about  0.85  MHz. 
Nevertheless,  the  value  of  h;)  obtained  herein  vields  values  of  h  at  the  higher  fre- 
quen  u  s  onsisti  r,t  w  ith  those  ol  Reference  2. 
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1'he  calculated  values  of  true  height  by  the  present  method  for  the  Boulder  and 
Garchv  ionograins  are  given  in  Figures  4  and  5.  The  calculated  values  of  h^  are 

the  initial  values  for  these  curves.  For  comparison,  the  previously  calculated 

■>  3 

values  of  true  height  for  Boulder”  and  for  Garchy  are  also  plotted  in  Figures  4 
and  a.  For  both  cases,  the  true  height  values  by  the  present  method  are  in  good 
agreement  with,  and  lying  slightly  below,  the  corresponding  curves  of  References  2 
and  3.  The  procedure  presented  herein  accordingly  appears  to  offer  a  simple 
method  without  undue  computational  labor  for  obtaining  true  height  within  the 
accuracy  of  the  experimental  data. 
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Figure  !.  Comparison  of  Real  Height  by  Present  Method  for  Boulder 
lonogram  itli  Previous  Calculation 
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Figure  5.  Comparison  of  Real  Height  by  Present 
Method  for  Garchy  Ionogram  With  Previous  Cal¬ 
culation 


5.  REPRESENTATION  OF  VIRTUAL  HEIGHT  DATA  IN  SIMPLE 
PARAMETRIC  FORM 

Although  the  calculations  of  the  increment  in  real  height  is  simple  in  principle, 
it  does  require  the  numerical  evaluation  of  an  integral  for  each  probing  frequency. 
We  have  therefore  explored  the  possibility  of  representing  the  virtual  height  by  a 
simple  parametric  model  from  which  the  real  height  can  be  rapidly  determined. 
Calculations  over  the  range  of  data  considered  here  showed  that  a  one-parameter 
model  was  sufficient  to  represent  the  data  with  good  accuracy.  However,  in  going 
from  the  lower  range  of  frequencies  to  the  higher  range,  it  was  found  necessary 
to  change  the  parametric  form  in  a  simple  way  in  order  to  maintain  the  same 
accuracy. 

Since  the  calculations  will  be  performed  for  several  values  of  the  probing  fre¬ 
quency  fy,  it  is  found  advantageous  to  utilize  the  normalized  probing  frequency 
x  =  f/fv-  Plots  of  the  increment  in  virtual  height  for  two  values  of  iy  are  presented 
in  Figures  6  and  7.  The  plots  have  been  normalized  to  the  maximum  increment 
in  virtual  height.  For  convenience  in  performing  the  calculations,  since  the  low 
frequencies  contribute  very  little  to  the  calculated  height,  the  curves  have  been 
extrapolated  to  zero  excess  height  at  zero  frequency.  It  may  be  noted  at  this  point 
that  the  higher  frequency  curves  show  considerably  more  curvature  than  those  at 
the  lower  frequencies.  It  is  this  increase  of  curvature  with  frequency  that  neces¬ 
sitates  the  change  in  parametric  form  in  going  to  the  higher  frequencies. 


Figure  6.  Plot  of  Virtual  Height  Function  g  Against 
Frequency  x  for  Boulder  Ionogram;  the  Function  g 
is  Normalized  to  Maximum  Virtual  Height 


Figure  7.  Plot  of  Virtual  Height  Function  g  Against 
Frequency  x  for  Garchy  Ionogram;  the  Function  g  is 
Normalized  to  Maximum  Virtual  Height 


In  looking  for  a  possible  parametric  fit  for  these  curves,  we  note  that  they  are 
approximately  linear  in  the  lower  range  and  then  bend  up  rapidly  as  the  frequency 
increases.  We  have,  therefore,  tried  a  polynomial  form  consisting  of  a  linear 
term  and  a  higher  order  term,  of  the  form 


h'  -  h' 


=  a ,x  + a„x 


where  h^^  is  the  maximum  value  of  the  virtual  height  and  h^  is  the  extrapolated 
value  of  h'  to  zero  frequency.  A  numerical  study  indicated  that  the  data  in  the 
lower  range  of  frequencies  could  be  fitted  by  the  polynomial 

g^  =  ax  +  (1  -  a)xl  (17) 

while  for  the  higher  range  of  frequencies  the  form 

g.,  =  bx  +  (1  -  b)x6  (18) 

was  a  good  model  for  the  data.  Plots  of  Eqs.  (17)  and  (18)  for  several  values  of 
the  governing  parameter  a  or  b  are  given  in  Figures  8  and  9.  To  show  the  kind  of 
fit  obtained,  several  data  points  from  Figures  6  and  7  are  replotted  on  these  fig¬ 
ures.  A  comparison  of  the  ionogram  data  and  the  model  curves  indicated  that  the 
former  could  be  represented  by  the  parametric  values: 
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Utilizing  Eq.  (2)  and  performing  the  required  integrations  we  obtain  for  the  real 
height  h, 

lower  frequencies,  Eq.  (17) 

h-h0“P<hmax‘h0,'“I0  +  (1-a,I1*  (li 

and  higher  frequencies,  Eq.  (18), 

h  ‘  ll0  =  |T(hmax  ’  h0)|bI0  +  (1  *  b)12>  *  (2‘ 


Figure  8.  Virtual  Height  Function  g  for  Lower  Frequen¬ 
cies  (Eq.  17)  as  a  Function  of  Frequency  for  Several  Val¬ 
ues  of  Parameter  a;  Data  Points  at  the  Lower  Frequen¬ 
cies  From  Both  lonograms  are  Also  Plotted  on  the  Figure 


Figure  P.  Virtual  Height  Function  g  for  Higher  Frequen¬ 
cies  (Eq.  18)  as  a  Function  of  Frequency  for  Several  Val¬ 
ues  of  Parameter  b;  Data  Points  at  the  Higher  Frequen¬ 
cies  From  Both  lonograms  are  Also  Plotted  on  the  Figure 


where  p  is  an  u>’t a'a^:  value  of  p.  Tno  i  integrals  arc  given  ov 


Evaluation  oi'  these  equations  yielded  the  real  heights: 
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As  seen  from  f  igures  5  and  0,  these  values  are  in  good  agreement  with  those  ob¬ 
tained  by  numerical  integration.  For  the  lower-frequency  Boulder  case.  Figure  8, 
the  choice  oi  a  proper  value  of  the  parameter  a  is  not  as  clear  cut  as  for  the  other 
cases.  However,  this  need  not  concern  us  greatly  since,  as  shown  by  an  incre¬ 
mental  analysis,  a  change  of  0.  On  in  the  parameter  a  results  in  a  change  of  only 
about  1  km  in  the  real  height.  The  parametric  procedure  presented  here  thus 
appears  to  be  a  useful  alternative  to  numerical  integration  in  tin  determination 
of  real  height  by  the  present  method. 

6.  SI  MM  AIM  AMI  CONCLUDING  REMARKS 

A  direct  method  for  determining  real  height  from  virtual  height  data  has  been 
applied  to  two  ionograms,  corresponding  to  data  obtained  at  Boulder,  Colorado, 
and  Garchy,  France.  The  initial  height  was  obtained  by  adapting  the  present 
method  to  the  one-start  procedure  developed  by  Howe  and  McKinnis  for  this 
problem.  The  results  obtained  for  the  two  ionograms  were  in  good  agreement 
with  those  obtained  from  previously-reported  accurate  numerical  calculations. 
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to  obviate  the  numerieal  integration  required  and  thus  make  the  present 
mi  tiled  more  rapid  and  convenient,  a  procedure  based  on  the  modeling  ot  the 
virtual  height  uata  in  terms  ot  a  polynomial  governed  by  a  single  parameter  lias 
been  developed.  The  parameter  may  be  determined  by  fitting  the  model  to  the 
given  data.  The  results  obtained  with  the  model  were  in  good  agreement  with 
those  previously  obtained  by  direct  use  of  the  data. 

The  present  method  affords  a  simple  and  convenient  means  of  obtaining  real 
heights  from  virtual  height  data,  within  the  accuracy  of  the  experimental  data. 


